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We investigate two classes of models of quintessential inflation, based upon canonical as well as 
noncanonical scalar fields. In particular, introducing potentials steeper than the standard expo¬ 
nential, we construct models that can give rise to a successful inflationary phase, with signatures 
consistent with Planck 2015 results. Additionally, using nonminimal coupling of the scalar field with 
massive neutrino matter, we obtain the standard thermal history of the Universe, with late-time 
cosmic acceleration as the last stage of evolution. In both cases, inflation and late-time acceleration 
are connected by a tracker solution. 
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I. INTRODUCTION 

Theoretical and observational consistency demands 
that the standard model of the Universe should be com¬ 
plemented by an early phase of rapid expansion dubbed 
inflation [1-14], as well as by late time cosmic accel¬ 
eration [15-20]. The latter is now accepted as a phe¬ 
nomenon of nature supported by independent sets of ob¬ 
servations, whereas inflation still awaits similar confirma¬ 
tion. The relic gravitational waves generated quantum- 
mechanically during inflation would have been a clear and 
direct signal of inflation [21-25]. In case the large value 
of r, investigations of B-mode polarization [26] could be¬ 
come a powerful tool to falsify the inflationary paradigm. 
Unfortunately, the Planck 2015 results [27, 28] seem to 
further shrink the bound on the tensor-to-scalar ratio of 
perturbations such that r = 0 is not ruled out [27-29]. 

Needless to say that inflation is one of the most beauti¬ 
ful and simple idea that not only resolves the inconsisten¬ 
cies of the hot big-bang such as the flatness problem, the 
horizon problem and others, but also provides us with a 
mechanism of generation of primordial perturbations. As 
for late time cosmology, the standard model of the Uni¬ 
verse is faced with yet another problem related to the age 
of the Universe, which is a late-time phenomenon [30-32], 
Interestingly, the resolution of the inconsistency within 
the framework of standard lore, asks for late-time cosmic 
acceleration, which was indeed confirmed directly by la 
supernovae observations in 1998 [33, 34] and was indi- 
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rectly supported by other probes independently [35, 36]. 
Obviously, accelerated expansion plays an important role 
in the history of Universe, both at its early and late 
stages. 

Often, these two regimes of accelerated expansion are 
treated independently. However, it is both tempting and 
economical to think that there is a unique cause respon¬ 
sible for both phases of acceleration a la quintessential 
inflation [37-64], which refers to unification of both con¬ 
cepts using a single scalar field. Consistency of the sce¬ 
nario demands that the new degree of freedom, namely 
the scalar field, should not interfere with the thermal 
history of the Universe, and thereby it should be “invis¬ 
ible” for the entire evolution and reappear only around 
the present epoch giving rise to late-time cosmic acceler¬ 
ation. It is, indeed, challenging to build a model which 
could successfully comply with the said requirements. 

First of all, one needs to construct an inflationary 
phase with a successful exit. Furthermore, in this sce¬ 
nario one needs an alternative reheating method (since 
the scalar field must survive till late times the conven¬ 
tional reheating is not applicable) and instant preheat¬ 
ing [65-67] is one of the efficient mechanisms that allows 
conversion of a part of the scalar field energy into radia¬ 
tion. In the post inflationary era till the present epoch, 
the field potential should be steep, allowing the radiation 
domination to commence, followed by a thermal history 
as envisaged by hot big bang. The latter is necessary for 
sending the field into hiding after the end of inflation. 
In particular, the post inflationary dynamics is charac¬ 
terized by a field that evolves into the kinetic regime for 
quite some time, but it then overshoots the background 
and gets frozen on its potential due to Hubble damp¬ 
ing. As the background energy density redshifts to the 
order of the field energy density, the field resumes its 
evolution. In case the potential is of a steep exponential 
form or steeper, the field tracks the background until late 
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times [64]. In case the potential is effectively shallow at 
late times, the field would exit from the scaling regime 
to slow roll. 

These features look very viable and pleasing, since it 
is implied that the late time evolution is broadly inde¬ 
pendent of initial conditions. The main reason for de¬ 
manding tracker behavior [68] after inflation is related to 
the hope of alleviation of the fine tuning. However, if we 
consider the interaction of the scalar held with matter, 
the mass of the scalar is destabilized, bringing back the 
same level of fine tuning with the cosmological constant 
paradigm [64]. 

In this paper we shall investigate models of quintessen¬ 
tial inflation using canonical (Sec. II) as well as non- 
canonical fields (Sec. Ill) with tracking behavior. In par¬ 
ticular, we are interested in constructing models that can 
produce a successful inhationary phase (SubSec. IIA for 
canonical held and SubSec. Ill A for noncanonical held), 
with signatures consistent with the Planck 2015 results, 
and then lead to the standard thermal history of the Uni¬ 
verse, with late-time acceleration as the last stage (Sub¬ 
Sec. IIB for canonical held and SubSec. IIIB for non¬ 
canonical held). Finally, in Sec IV we summarize our 
results. 


II. UNIFYING INFLATION AND 
QUINTESSENCE USING A CANONICAL 
SCALAR FIELD 


In this section we study quintessential inflation using 
a canonical scalar held. We consider the action 


S = J d 4 Xy/^g 

+5 m + <S r , 


M 2 1 


(1) 


with Mp\ the Planck mass, <f> the scalar held, and V((f>) 
its potential. In the above action we have additionally 
considered the matter and radiation sectors S m and S r 
respectively. These sectors can be neglected at the in¬ 
hationary stage, however they will gradually play an im¬ 
portant role, giving rise to the standard thermal history 
of the Universe and finally to the late-time accelerating 
phase. As usual, we focus on the case of a hat Friedmann- 
Robertson-Walker (FRW) geometry, with metric 

ds 2 = — dt 2 + a(t) 2 5ijdx l dx :i , (2) 

where a(t) is the scale factor. Friedman equations are 
given by 

3 H 2 Mpi = p m + p r + — (j ) 2 + V ((/)) 

( 2 H + 3 H 2 ) M 2 x = -± Pl - 2 + , ( 3 ) 

and the equation of motion for the scalar held has the 
standard form 

dV 

(f> + 3H<f> + —— = 0 . (4) 


A. Inflation 


In what follows, we shall first analyze the inflationary 
phase in this scenario, focusing on the signatures on the 
observables that allow for a comparison with the Planck 
data. As usual, in the inflationary phase one may neglect 
S m , <S r and <S„, and thus the dynamics of inflation, as well 
as its observational signatures, are determined solely by 
the scalar held and its potential. In particular, given the 
potential one introduces the slow-roll parameters 
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Additionally, the usual condition for ending inflation is 
simply 


t* —0end 


= 1, 


( 8 ) 


where the subscript end represents the value at the end of 
inflation (we follow the same convention in the rest of the 
paper). The number of e-foldings is calculated through 
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(9) 


Hence, observables like the tensor-to-scalar ratio (r), 
the scalar spectral index (n s ) and its running (a s = 
dn s /dlnfc), can be written as 


r ~ 16e, (10) 

n s « 1 — 6e + 2 ? 7 , (11) 

a s « 16e?7 — 24e 2 — 2£ 2 . (12) 


Keeping in mind the discussion in the introduction, we 
consider the potential 

V = V 0 e~ xr/M ^, (13) 

where Vo and A are the usual parameters. Note that com¬ 
pared to standard exponential potential, we have allowed 
for one more parameter, namely n, which would influence 
the steepness of the potential. The case n = 1 has been 
extensively studied in the literature [25, 40, 61, 69-71] 
and thus in the following we consider the case n ^ 1. 
Moreover, we consider the cases n / 2 and n = 2 sepa¬ 
rately, since the corresponding expressions are different 
in these cases. 


1. n ± 2 


In this case the slow-roll parameters (5)-(7) have the 
following form, 
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rj = -M^ 2n nX(p n - 2 [Mpj(n - 1) - n\cj> n } , (15) 


$ 2 = Mpj - 4 n n 2 A 2 </> 2n-4 [M|f(n 2 - 3n + 2) 

-3Mft(n - l)n\cj) n + n 2 A 2 ^> 2 "] , (16) 

where <j> is the value of the field at the horizon crossing. 
Additionally, condition ( 8 ) gives 


4>end = Mpi 


n 2 X 2 


(17) 


and thus from (9) we obtain 

M = 


M n ~ 2 

P1 - «£r) 


nA(n — 2) 
1 

nA(n — 2) 


V _T * - 


n 2 A 2 


1 — n ~| 
2n — 2 


.(18) 


in order to obtain agreement with Planck Collaboration 
[27, 28], 

In Fig. 1 we depict the allowed regions in the n — X 
parameter space that can give n s = 0.9644 ± 0.0049 [28] 
and 0 < r < 0.149 for M = 60. We clearly see that the 
parameter n must be larger than 5. It is interesting to 
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One can revert this expression in order to get 
$ = Mpi Q(n, A, A/”) 


(19) 


FIG. 1: The shaded region marks the allowed region in the 
n — X parameter space that can lead to n s = 0.9644 ± 0.0049 
and 0 < r < 0.149 for Af = 60. 


with 


Q(n,X,JV) = jnA 

( 20 ) 

which allows to eliminate <j) in favor of Af in the slow-roll 
parameters. In particular we acquire: 

e=in 2 A 2 Q(n,A,A0 2 "- 2 , (21) 


notice that if we exclude the zero value, for instance if 
we consider r > 0 . 01 , then the corresponding region is 
significantly reduced. In Fig. 2 we depict the regions in 
the n — X parameter space that can give n s = 0.9644 ± 
0.0049 and 0.01 < r < 0.149 for Af = 60 and A f = 70, 
where the aforementioned feature is clear for A f = 60. 
However, this does not seem to be the case according 
to both 2013 [72] and 2015 [27, 28] Planck data sets. 
Fig. 2 also shows that the parameter space increases if 
we increase the value of A7 from 60 to 70. 


(n — 2)Af + nA2 5 F^ r T) 


2 A 2 


2 (" — 1 ) 


77 = nXQ(n, A, A/ - )” 2 {1 — n + nXQ(n, X,Af) n } , ( 22 ) 

£ 2 = n 2 A 2 Q(n, A, A0 2n ~ 4 

• {2 + n [n — 3 — 3 (n — 1)A Q(n, A, A f) n 

+nX 2 Q(n, X,J\f) 2n ]} . (23) 

Thus, the tensor-to-scalar ratio (r), the scalar spectral 
index (n a ) and its running a s , can be calculated straight¬ 
forwardly as functions of n, A, AT using ( 10 )-( 12 ). 

Let us now use the above expressions to determine for 
which combinations of n, A and e-folding Af we obtain 
values of n s and r in agreement with the Planck 2015 
results. In particular, we desire to obtain n s = 0.9644 ± 
0.0049 (68 % confidence level, Planck TT,TE,EE+lowP) 
consistent with the Planck 2015 results [28] and 0 < r < 
0.149 (recent joint analysis of BICEP2/Keck Array and 
Planck data gives r 0 .05 < 0.12 at 95% confidence [29] 
and when running of the scalar spectral index is allowed 
Planck 2015 results give r < 0.149 [28] at 95 % confi¬ 
dence). As a starting point, and for completeness, we are 
interested in obtaining 0 < r < 0.149, and thus describ¬ 
ing the limiting cases of both Collaborations. However, 
later on we will focus on the low values of this range, 
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FIG. 2: The blue shaded region (upper shaded region) and the 
green shaded region (lower shaded region) mark the allowed 
region in the n — X parameter space that can lead to n s = 
0.9644 ±0.0049 and 0.01 < r < 0.149 for M = 60 and Af = 70 
respectively. 

In order to investigate further the effect of the param¬ 
eters n and A, for different e-folding number A7, on r and 
n s , we include various figures. Firstly, in Fig. 3 we depict 
r versus A, for different values of n and e-folding AT, while 
in Fig. 4 we show r versus n, for fixed A and different Af. 
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FIG. 3: The tensor-to-scalar ratio r versus the parameter A, 
for different values of the parameter n and e-folding Af. 


FIG. 6: The scalar spectral index n B versus the parameter n, 
for fixed parameter A and different e-folding value Af. 
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FIG. 4: The tensor-to-scalar ratio r versus the parameter n, 
for fixed parameter A and different e-folding value Af. 
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FIG. 7: la (yellow) and 2a (light yellow) contours for Planck 
2015 results (TT + lowP + lensing + BAO + JLA + Ho) [27], 
and la (grey) and 2a (light grey) contours for Planck 2013 
results (Planck + WP+BAO) [72] (note that the la region of 
Planck 2013 results is behind the Planck 2015 results, hence 
we mark its boundary by a dotted curve), on n B — r plane. 
Additionally, we depict the predictions of our scenario, for 
varying A (between IIP 6 and 10 -3 ,), and n being 4 or 6, with 
the e-folding value Af being 50 or 70. 


FIG. 5: The scalar spectral index n s versus the parameter A, 
for different values of the parameter n and e-folding Af. 


Similarly, in Fig. 5 we depict n s versus A for different 
values of n and Af, while in Fig. 6, we show r versus n, 
for fixed A and different Af. 

It is clear from the above discussion that the scenario 
at hand, with the potential (13), can give rise to n s and r 
in agreement with both the Planck 2013 results [72] and 
the Planck 2015 results [27, 28]. In order to present these 
features in a more transparent way, in Fig. 7 we depict 
the predictions of our scenario for varying A, and n being 
4 or 6, with the e-folding value Af being 50 or 70, on top 


of the la and 2<r contours of the Planck 2013 results [72] 
as well as of the Planck 2015 results [27]. As we observe, 
as n or Af increase, the predictions move towards the 
core of the data. Furthermore, in Fig. 8 we present the 
corresponding situation, but for varying n, and A being 
10~ 4 or 10“ 5 , with the e-folding value Af being 50 or 
70. As we observe, as Af increases the predictions move 
towards the core of the data. Hence, we deduce that the 
larger parametric freedom that was introduced by the use 
of the additional “steepness” parameter n, comparing to 
models with only the parameter A, can lead to the desired 
r-n s behavior. 

For completeness, let us make a comment on the pre- 
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FIG. 8: la (yellow) and 2a (light yellow) contours for Planck 
2015 results (TT + lowP + lensing + BAO + JLA + Ho) [27], 
and la (grey) and 2a (light grey) contours for Planck 2013 
results (Planck+WP+BAO) [72] (note that the la region of 
Planck 2013 results is behind the Planck 2015 results, hence 
we mark its boundary by a dotted curve), on n s — r plane. 
Additionally, we depict the predictions of our scenario, for 
varying n (between 4 and 20), and A being 10~ 4 or 1CP 5 , with 
the e-folding value Af being 50 or 70. 


diction of the scenario at hand on the running spectral 
index a s = dn s /d In k ~ 16ep — 24e 2 — 2£ 2 . Using (12) 
and (21)-(23), we can calculate it for various values of 
A, n and A f, and we present the results on the a s — n s 
plane in Fig. 9. On the same graph we depict the lcr and 
2er contours of the Planck 2013 results [72] as well as of 
the Planck 2015 results [27, 28]. As we observe, as n or 
A f increase the predictions move towards the core of the 
data, and especially for the parameter values of Figs. 7 
and 8 we obtain a remarkable agreement with the Planck 
2015 results [27]. 

Let us now calculate the energy scale of inflation us¬ 
ing the COBE normalized value of density perturbations, 
which can be represented by the following fitting function 
[73] 

S H (n s ,r) = 1.91 x 10" 5 e lol(1_ns) /Vl + 0.75r. (24) 


On the other hand, the scalar perturbation spectrum is 
given by 


A 2 s {k) 


V 

(150t r 2 M 4 l£ ) ’ 


(25) 


and at the horizon crossing (k = fc* = a*Ff*) it becomes 


A 2 (fc») = 7"»*“ 1 4 . (26) 


Using Eqs. (24), (25) and (26), we can have the estima¬ 
tion of some model parameters. For instance, for r = 0.05 


FIG. 9: la (yellow) and 2a (light yellow) contours for Planck 
2015 results (TT,TE, EE + lowP) [27], and la (grey) and 
2a (light grey) contours for Planck 2013 results (ACDM + 
running + tensors) [72], on a a — n s plane. Additionally, we 
depict the predictions of our scenario, for varying A (between 
10 -6 and lCF 3 ), n = 6, AT = 70, for varying n (between 4 
and 20), A = 1CF 4 , AT = 70, and for varying n (between 4 
and 20), A = 10 -4 , Af = 50. 


(best fit value of r according to [29]), J\f = 70 and n = 6, 
equations (10) and (21) give A = 1.46 x 10 -9 , which leads 
to Vo = 3.39 x 10~ 9 Mpp Additionally, for the same val¬ 
ues of r, n, A and Vo, the value of the potential at the 
commencement of inflation is Un = 1.4 x lO^Mpp which 
provides the scale of inflation as = 1.49 x 10 16 GeV. 

Finally, let us discuss on the constraints on reheating 
temperature from relic gravitational waves. As shown in 
Refs. [24, 25, 61], the ratio of the energy densities of relic 
gravitational waves produced during the kinetic regime 
(p g ) and radiation (p r ), is 


P± 

Pr 


end 


3?r (Ps\ 

64/iow \pJ 


(27) 


where “eq” represents the equality of radiation and scalar 
field energy densities. Moreover, the square of relic grav¬ 
itational wave amplitude writes as 


h 


2 _ 
GW — 


Hi 


87rMp[ ’ 


(28) 


where H ln is the Hubble parameter at the commencement 
of inflation, which is « ^Vi n /(3Mp\). 

Nucleosynthesis imposes a constraint on the ratio of 
the relic gravitational waves and radiation energy densi¬ 
ties, namely (p g /p r )eq < 0.01 [36]. Hence, this provides 
the constraint on the amount of radiation energy density 
at the end of inflation, that is 


9U 0 2 -A[er(n,A,A0+(x|) n/(n 4) 

Pr,en d „ ^ e 


( 29 ) 
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Furthermore, the temperature at the end of inflation is 
Fgnd = p] {td• Therefore, using Eq. (29) we can also get a 
constraint on the temperature at the end of inflation, that 
is the reheating temperature. If we consider r = 0.05, 
AT = 70 and n = 6 , then we have already seen that 
A = 1.46 x 10 -9 and Vq = 3.39 x lO^Mpj. Hence, for 
these values of the model parameters we get 

r end > 2.264 x 10 14 GeV. (30) 

2. n = 2 

In this paragraph we present the results in the n = 2 
case for completeness. In this case, the slow-roll param¬ 
eters (5)-(7) become 


p = 2e — 2A 

£ 2 = 4e (e — 3A) . (31) 

Furthermore, the number of e-foldings is 

"-**{*?)• (32) 

with </>end = ^ 2 \- Thus, <f> can be expressed through 
A and N as <f> = ^^ e_2AA ^i and therefore we can write 

e = e _4AA ^. Hence, the tensor-to-scalar ratio r, the scalar 
spectral index n s and its running a s , are written as 

r « 16e = 16e~ 4AM , (33) 

n s « 1 - 6 e + 2r/ = 1 - 2e~ 4XM - 4A, (34) 

A s « 16er 7 - 24e 2 - 2 £ 2 = - 8 Ae~ 4AAA . (35) 

Unfortunately, as one can see expressions (33), (34), 
(35) cannot lead to values in agreement with Planck re¬ 
sults for 50 < M < 70, independently of the A value. 
Thus, we do not investigate this case in more detail. 

B. Late Time Dynamics 

In this subsection we investigate the late-time behavior 
of the above scenario. For usual steep exponential poten¬ 
tial (n = 1 ) we know that during the post-inflationary 
dynamics the scalar field rolls down the potential, and 
its energy density scales as p$ ~ a” 6 . Due to the in¬ 
creased Hubble damping, the scalar field stops evolving, 
so its energy density eventually becomes comparable to 
the background, and it again starts evolving and scales 
with the background up to late times, thus leaving no 
place for late-time acceleration. This class of solutions is 
known as scaling solutions [74]. In Fig. 10(a) we present 
such a scaling behavior of the scalar field energy density 
for an exponential potential. 
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(b) 

FIG. 10: Evolution of the energy densities of matter (dotted 
green), radiation (dashed blue) and scalar field (solid red), as a 
function of the redshift (z = ao/a — 1 with oq = 1 the present 
scale factor), in the case of the minimally coupled scenario 
(1). The upper graph is for n = 1 and A = 4, while the lower 
graph is for n = 2.5 and A = 4. 


Now let us see what happens when the potential is 
steeper than the exponential one, i.e. the case where 
n > 1. Similarly to the exponential potential, in this 
case too the energy density of the scalar field decreases 
rapidly, and due to the large Hubble damping the scalar 
field stops evolving and thus eventually its energy den¬ 
sity becomes comparable with the background one. But 
unlike the exponential case, now, due to the very steep 
nature of the potential, the scalar field cannot follow 
the background during the high redshift and thus again 
Pcf, ~ a~ e . This results in a rapid decrease in the scalar 
field energy density and hence again the scalar field ex¬ 
periences large Hubble damping due to the background, 
and therefore it repeats the same behavior as explained 
earlier (see Fig. 10(b)). However, we mention that this 
can happen only during large redshifts, where the field 
value is not so large. 

On the other hand, at late times, when the field evolves 
to a large value, the picture is different. In order to 
understand the behavior during late times, let us consider 
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the function F = V"V/V' 2 . For an exponential potential Additionally, the continuity equation for massive neu- 
r = 1, however for the steeper potential (13) we have, trinos is given by 


F = 1 - 


(n - 1) Mg t 
nX <f> n 


(36) 


Pv + 3 H (p v + p v ) = 7 (p v - 3 p v ) —— 

M-pi 


(43) 


From this expression it is clear that for large <f> and n > 1 
the function F approaches 1, i.e., for asymptotically large 
field values, the nature of the potential eventually be¬ 
comes similar to the exponential one. Thus, if at late 
times the field value is sufficiently large, we obtain a scal¬ 
ing solution, and indeed Fig. 10(b) confirms this. 

Unfortunately, as can be also seen in Fig. 10(b), we 
cannot obtain late-time acceleration for the potential 
(13). In order to achieve late-time acceleration we need 
a mechanism to exit from the scaling behavior, that is 
to obtain the tracker behavior [68] at late times. For 
this purpose we can consider a nonminimal coupling 
between the scalar field and massive neutrinos, as in 
Refs. [61, 62, 75] (also see Refs. [76-90]), and we start 
with the action: 


S= d 4 x\/^g 


M 2 j „ 1 


2 2 
+ iS r + iSj, ( C 2 g a p , H/,,) , 




(37) 


where 


C 2 _ ^2 e 270/Mpi . 


(38) 


In flat Friedmann-Robertson-Walker (FRW) cosmology 
the two Friedmann Eqs. (3) modify to: 

3iF 2 Mpj = —(j) 2 + V (</>) + p m + pr + Pu (39) 

(2 H + 3 H 2 ') M 2 x = ~\tf + V(4>) - l -Pr - Pv, (40) 
while scalar-field equation (4) now becomes 

<41) 

Before proceeding further we should mention here that 
massive neutrinos are relativistic (p v = p v / 3) for the 
most of the expansion history of the Universe and become 
nonrelativistic (p„ = 0) only after the redshift znr € 
(2 — 10) for neutrino mass range m u £ (0.015 — 2.3) eV 
[76, 77]. So concerning the neutrino equation-of-state 
parameter we shall consider the following ansatz [61]: 


w v {z) 


Pv_ 

Pv 


1 

6 


1 + tanh 


ln(l + z) - z eq 
Zdur 


(42) 


where z eq and Zdur are two parameters which determine 
the redshift around which the transition of w u from 1/3 to 
0 starts and how fast the transition happens, respectively 
and the values of these two parameters depend on the 
redshift znr- 


The last term of Eqs. (41) and (43) is effectively zero 
when neutrinos behave like radiation, however it becomes 
non-zero when neutrinos become nonrelativistic (p v = 
0). So the nonminimal coupling between the massive 
neutrinos and the scalar field affects the expansion of 
the Universe only during the late times and as can be 
deduced from (41), an effective potential forms, which 
reads as 

U eff = V(0) + pp O e T(0 - 0o)/Mpi , (44) 

where (f > 0 and p v 0 are the present values of the field 
and of the massive neutrino energy density, and p u = 
Pv. This effective potential clearly has a 
minimum for 7 > 0, which forms at late times. Hence, 
the scalar field oscillates around this minimum, and even¬ 
tually it settles down to the minimum as the oscilla¬ 
tions decrease due to the Hubble friction. For clarity, 
in Fig. 11(a) we depict the numerically evolved effec¬ 
tive potential, normalized by the present critical density 
(p c o), around its minimum. It should be noted that the 
minimum value of the effective potential I 4 ff,min normal¬ 
ized by the present critical density p c 0 is ~ 1 , which 
implies that pde ~ I4ff jln i n since p c 0 « Pde- Moreover, 
from Fig. 11(b) we can see that T4ff,min ~ Pv, min, where 
Pv,min is the energy density of the massive neutrinos at 
the minimum of the effective potential. Furthermore, 
from Fig. 11(c) we observe that p„. m i n ~ p v o- Hence, in 
summary we can deduce that in the model under con¬ 
sideration the dark energy scale is related to the present 
energy density of the massive neutrinos. 

As we described above, the nonminimal coupling be¬ 
tween the scalar field and the neutrinos and the induced 
effective potential, is adequate to lead to late time accel¬ 
eration. Indeed, in Fig. 12 we depict the evolution of 
the various energy densities, and we can clearly see the 
tracker behavior of the scalar field and the onset of the 
dark-energy dominating phase. 

In order to present the thermal history of the Universe 
in a more transparent way, we introduce the dimension¬ 
less density parameters for matter, radiation, neutrinos 
and scalar field, respectively given by 

o - _ Em _ 

3 H 2 M 2 X ’ 

ci - _ El _ 

r ~ 3H 2 M 2 ! ’ 

ci - _ El _ 

3H 2 Mp l ’ 

O — ^ 

0 3H 2 Mp l ’ 

where p$ = (1/2 )<p 2 + V, and in Fig. 13 we depict their 
evolution as a function of the redshift. As we observe, 


(45) 

(46) 

(47) 

(48) 
















54.54 54.55 54.56 54.57 54.58 54.59 54.60 

<P 


(a) 



(b) 



(c) 

FIG. 11: Top: The minimum of the effective potential ()f) for 
pvo/pco = 0.0054 (Planck 2015 results gives Tl v oh 2 < 0.0025 
and Ho = 67.74 ± 0.46 km s -1 Mpc -1 [27]), 7 = 800 and 
A = 10 -8 , with p c o the present critical density. Middle: The 
ratio of the minimum of the effective potential over the mas¬ 
sive neutrino energy density versus the field value around the 
minimum of the effective potential. Bottom: The neutrino 
energy density, normalized with its present value, versus the 
field value around the minimum value of the effective poten¬ 
tial. For all the plots 7 = 800, A = 10 -8 , n = 6 , z eq = 
2.55 and Zdur = 3. 


we can reproduce the thermal history of the Universe, 
starting from a scalar field kinetic regime, then enter¬ 
ing the radiation and matter regimes, and finally result¬ 
ing to late-time dark-energy dominated era. Finally, for 
completeness, in Fig. 14 we depict the corresponding 
behavior of the scalar-field equation-of-state parameter 
Wcf, = P4,/ Pcj,, as well as of the effective (total) equation- 



FIG. 12: Evolution of the energy densities of matter (dotted 
green), radiation (blue short dashed), scalar field (red solid) 
and massive neutrinos (purple long dashed), as a function of 
the redshift, in the case of the nonminimally coupled scenario 
(37), for 7 = 800, A = 10 “ 8 , n = 6 , z eq = 2.55 and z dul = 3. 

of-state parameter w e s = Ptot/ptot = —1 — 2H/3H 2 . 
From this figure we verify that around the present era 
the potential term dominates over the kinetic one in the 
scalar-field energy density, which leads to be around 

- 1 . 



Log,,, (1 + z) 


FIG. 13: Evolution of the density parameters of radiation 
(blue dashed), matter (green dashed-dotted) , scalar field (red 
dotted) and neutrinos (black solid), respectively, in the case 
of the nonminimally coupled scenario (37), for 7 = 800, A = 
10 ” 8 , n = 6 , z eq = 2.55 and z d ui = 3. 


III. UNIFYING INFLATION AND 
QUINTESSENCE USING A NONCANONICAL 
SCALAR FIELD 

In this section we shall be interested in constructing 
quintessential inflation using a class of models with non- 
canonical scalar field [61, 64, 75, 91 94] (see also Refs. 
[95-98] for unification using noncanonical phantom field). 
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FIG. 14: Evolution of the scalar-field, (blue doted) and effective 
(total) (red solid) equation- of-state parameters, in the case of 
the nonminimally coupled scenario (37), for 7 = 800, A = 
1CP 8 , n = 6 , z eq = 2.55 and Zdur = 3. 


by V(a) ~ g-atr/Mpi, an( j thereby a should be small to 
comply with inflation [62], Concerning the parameter 
/3, it can be fixed by COBE normalization [62]. Hence, 
once the post inflationary behavior is guaranteed by the 
specific form of the kinetic function, inflationary require¬ 
ments can be obtained through appropriate choices of a. 

It is clear from the aforesaid that (49) can give rise to a 
viable model of quintessential inflation by adding a non- 
minimal coupling between the scalar field and the neu¬ 
trinos described in subsection IIB. In the following sub¬ 
section we shall describe inflation using the noncanonical 
model (49), and we will derive constraints on the model 
parameters in the light of Planck 2015 results [27, 28]. 
Finally, in a separate subsection we will examine the late¬ 
time, post-inflationary evolution. 

A. inflation 


In this case, naturally, we have tracking behaviour in 
the post inflationary era [61]. Comparing with thawing, 
the tracking behavior imposes tough restrictions on the 
post-inflationary evolution of the scalar-field dynamics, 
namely the field should mimic the background for most 
of the Universe history, and only at late times it should 
exit to slow roll regime. The latter is realized only for spe¬ 
cific potential forms, otherwise the field exhibits thawing 
behavior. In summary, in general it is difficult to acquire 
the tracking features after inflation, using simple poten¬ 
tials. However, the picture changes if we use a scalar 
field with noncanonical kinetic energy, since in this case 
it is easy to control the post-inflationary dynamics in the 
desired way. 

Let us consider the following action [61, 75]: 


5 = 


M ^ R _ _ V (fi) 


j d 4 xy^ 

+S m + S r + Sv(C z (4>)g a p', Tp), 


2 2 
-ill 


(49) 


with 


According to the above discussion, the noncanonical 
field action (49) can lead to inflation which commences 
in the small-field region. In particular, the slow-roll pa¬ 
rameters can be expressed as [62] 


e = 


Mp 1 


IdUV 


V 

e 


2 k 2 {(/>) \v d(/> J 

Mp\ de(4>) 


2 k 2 (<j)) 
-2 


= 2 e- 


= 2eq — 


a d<f> 
aMp\ dr/ 
k 2 d (j) 


~ e 


a 

Y 


~ 2ore, 


~y(l + J0,(53) 

(54) 

(55) 


with X = ffi 1 e a 't > l Mpl , and where we have used the ap¬ 
proximations a 1 and aCl, which should hold in the 
scenario under consideration. It is clear from (54) that 
inflation ends in the region with X 1, which quantifies 
the large-field approximation. In this approximation the 
number of e-foldings becomes [62] 


-^Inh+X- 1 ) , 
or v 7 


(56) 


k ~^ = ( a 2 ) 1 + /3 2 e^/ M Pi + 1 (50) 
V(<j>) = Mp ie _ “^ Mpi (51) 

C((j)f = Ce 2 ^ /Mpi , (52) 

where a, a, 7 and /? are the model parameters. The 
kinetic function k 2 (</>) has been suitably chosen according 
to the tracking post-inflationary requirements. 

In order to obtain the realistic intervals of the param¬ 
eter space, we start by noting that k(<j>) —> 1 in the 
large-field limit. Thus, in this case we obtain a canoni¬ 
cal field with exponential potential, whose slope is given 
by a, which should be large in order to adhere to nu¬ 
cleosynthesis constraint, namely a > 20 [62]. Similarly, 
in the small-field limit we can introduce a canonical field 
a = ( a/a)<f >, such that the field potential is approximated 


and thus it can be related to e through [62] 


e(AT) 


a 2 1 
2 1 - e -“ w ' 


(57) 


Let us note that the transition between small and large 
field regimes takes place when d 2 « 1 /J\f. Since inflation 
always ends in the region of large field, its commence¬ 
ment depends upon the range of inflation, which is in 
turn uniquely specified by the tensor-to-scalar ratio. In 
particular, a large value of r, or weak slow roll, would 
imply large field excursion. In that case inflation should 
commence around the boundary of transition, otherwise 
the commencement would be shifted to the large-field re¬ 
gion. 

The general expressions for r, n s and the running of 
spectral index a s , valid from small to large field regimes 
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n 

s 


In general the effective Lagrangian, after we fix the ig¬ 
norance, contains both renormalizable as well as non- 
renormalizable parts. As for renormalizable part, it in¬ 
cludes one loop corrections to classical Lagrangian. In 
our case, the latter is absent. The non-renormalizable 
part includes correction that are suppressed by inverse 
powers of the cut off. As long as we work quite below 
the cut off, we can safely use the classical framework. 
The effective Lagrangian in our case has the form, 


£ — £ci + 

i -1 


Ci0 n 


^-4+m 


di(d4>) 2 (/) m 4 \ 

^ — 4 +ni J 


(61) 


where we imagined series expansion of the potential and 
Ci,di are constants. In this case, marginal and relevant 
operators are absent. Thus, in case we work well below 
the cut off, we can safely ignore the correction and keep 
using the classical framework. 


FIG. 15: la (yellow) and 2a (light yellow) contours for 

Planck 2015 results (TT+lowP+lensing+BAO+JLA+Ho) 
[27, 28], and la (grey) and 2a (light grey) contours for Planck 
2013 results (Planck + WP+BAO) [36] (note that the la re¬ 
gion of Planck 2013 results is behind the Planck 2015 results, 
hence we mark its boundary by a dotted curve), on n s — r 
plane. Additionally, we depict the predictions of our scenario 
given by (58),(59), fora —> 0 and e-foldingAf varying between 
55 and 70. 


are given by 


r(Af, a) = 16e(A0 


8 a 2 


^ _g— ct 2 N 5 

n s (Af, a) = 1 — 6e + 2rj ss 1 — a 2 coth 
dn. 


m 

a, 4 


d In k 


= 16er] - 24e 2 - 2 £ 2 « - 


2 sinh 2 (^L) 


(58) 

(59) 
.(60) 


Hence, let us use these expressions in order to compare 
the predictions of the scenario at hand with the 2013 and 
Planck 2015 results. In Fig. 15 we depict the predictions 
of our model for a —> 0 and e-folding M varying between 
55 and 70, on top of the lcr and 2cr contours of the Planck 
2013 results [36] as well as of the Planck 2015 results 
[27, 28]. As we observe, the point for J\f = 55 lies outside 
the 2a contour of both Planck data sets, but higher val¬ 
ues of A( lie within the 2 a contour of both Planck data. 
Unfortunately, all predictions still lie outside the lcr con¬ 
tour of both data sets, and this becomes worse for larger 
values of a. 

As our potential (13) is of unusual form, a comment 
about the viability of model under quantum correction 
is in order. Indeed, since, V ~ e^ ,n > 5, the model 
would involve operator of dimensions higher than four 
if we imagine the series expansion of the potential. It 
then obviously raises the question whether the model 
would make sense if quantum corrections are invoked. 


B. Late Time Dynamics 

In the noncanonical scalar-field scenario at hand, after 
the end of inflation the Universe enters into a kinetic- 
energy-dominated regime, known as “kinetic regime”, 
and then it subsequently enters into the radiation, mat¬ 
ter and dark-energy eras [61]. The nonminimal coupling 
between the scalar field and massive neutrinos plays the 
main roll for the onset of late-time cosmic acceleration, as 
we analyzed in subsection IIB. In particular, when the 
massive neutrinos become nonrelativistic at late times, 
they contribute to the formation of the effective poten¬ 
tial with a minimum. As a result, the scalar field set¬ 
tles down to the minimum of the effective potential after 
damping of oscillations which ultimately gives rise to the 
late time acceleration. For a detailed dynamical analysis 
one can see Ref. [61]. The late-time attractor solution 
corresponds to an effective (total) equation-of-state pa¬ 
rameter given by [61] 


ffleff = - 


7 

1+7 


and to a scalar-field equation of state 


a 2 q(l + 7) 

3 + ct 2 7 (l + 7 ) ’ 


(62) 


(63) 


where a is the canonical scalar field which can be repre¬ 
sented in terms of the noncanonical scalar field <j> using 
the transformation [61] 

a = k((/>), (64) 

= (S) ■ <65) 


From (63) we can see that +7 = 0 , i.e. without a coupling 
between the scalar field and massive neutrinos, w a = 0 , 
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which implies that the scalar field will exhibit scaling be¬ 
havior even during late times, and will continue to follow 
the background even in the future. Hence, we do verify 
what we discussed earlier, namely that in the absence of 
the nonminimal coupling we cannot acquire late-time ac¬ 
celeration. On the other hand, from (62) we deduce that 
in order to obtain a de Sitter or nearly de Sitter solution 
( w e ff « — 1) we require 7>1. Thus, a large nonminimal 
coupling is needed in order to acquire late-time acceler¬ 
ation. 

Finally, note that the value of the effective potential at 
the minimum is directly proportional to the present neu¬ 
trino energy density [61, 64], Therefore, the dark energy 
scale is related to the neutrino energy scale, similarly to 
the analysis of Section II. Hence, the nonminimal cou¬ 
pling between the scalar field and massive neutrinos not 
only provides the late-time acceleration, but it addition¬ 
ally fixes the energy scale of dark energy. 

IV. CONCLUSIONS 

In this work we investigated two distinct classes of 
quintessential inflation, namely models based on a canon¬ 
ical scalar-field and models based on a noncanonical 
scalar field, where both scenarios exhibit tracking behav¬ 
ior. In both cases we considered a nonminimal coupling 
between the scalar field and the neutrinos, which is re¬ 
quired in order to trigger the late-time cosmic accelera¬ 
tion. 

In the canonical case we considered a potential of the 
form V ~ e x< ^ / Mpl , with n > 1, which has the prop¬ 
erty of slow roll near the origin but it becomes steep 
away from it. Hence, at early times, i.e. while the field 
is around the origin, this scenario can give rise to in¬ 
flation. Indeed, as we demonstrated in Figs. 1 and 2, 
for a range of the model parameters we can obtain a 
required phase of inflation with the spectral index n s 
and the tensor-to-scalar ratio r in a very good agree¬ 
ment with the Planck 2015 results [27, 28] and the joint 
analysis of BICEP2/Keck Array and Planck data [29]. 
Indeed, in Fig. 7 and Fig. 8 we showed that the pre¬ 
dictions of the scenario fall well inside the ler likelihood 
contours of both Planck 2013 results [72] and Planck 
2015 results [27, 28]. Additionally, for the representa¬ 
tive case of the parameter choice n = 6, A = 1.5xl0 -9 
and e-foldings A f = 70, the obtained r = 0.05 provides 
the estimate Vo = 3.39 x 10 _9 Mp[ or equivalently the 
scale of inflation, V^ 4 = 1.49 x 10 16 GeV. Finally, 
using nucleosynthesis constraints we obtained the lower 
bound on the temperature at the end of inflation, namely 
T end ^ 2.264 x 10 14 GeV. 

After the end of inflation, the steep potential derives 
the scalar field into the kinetic regime. Consequently, the 
field overshoots the background and freezes due to Hub¬ 
ble damping. The evolution of the field resumes soon af¬ 
ter the background energy density becomes comparable 
to field energy density. In the usual case of exponential 


potential (n = 1) the field follows scaling behavior. In the 
case of n > 1, the potential is steeper than the standard 
exponential and thus the field is driven away from the 
scaling track, which increases the Hubble damping lead¬ 
ing to the freezing of the field once again. As the field 
comes out of the freezing regime, its energy density red- 
shifts faster than the background, and this feature brings 
back the Hubble damping and so on. The said behavior 
keeps repeating till the field acquires large values. In that 
case F —> 1 and the field enters the regime which is an 
attractor, see Fig. 10(b). The latter allows us to obtain 
the subsequent radiation and matter regimes. Thus, we 
conclude that the scaling solution is also an attractor in 
case the potential is steeper than a standard steep expo¬ 
nential. To the best of our knowledge, this feature was 
not noted earlier in the literature. In case of generic val¬ 
ues of n and A, the field rolls in the domain of large <f> in 
the post inflationary era in which case the system enters 
into the tracking regime after the Hubble freezing ends, 
see Fig. 12. 

Finally, the nonminimal coupling between the scalar 
field and the neutrinos induces an effective potential, 
which leads the scalar field to drive the late-time ac¬ 
celeration. The larger the nonminimal coupling is, the 
deeper the minimum of the effective potential is, in which 
the field is settled after damped oscillations, exhibiting 
an equation-of-state parameter around —1 in the present 
epoch. 

As for the noncanonical scalar field, although both in¬ 
flation and the subsequent thermal history of the Uni¬ 
verse, including late-time acceleration, can be obtained, 
the specific values of the spectral index and of the tensor- 
to-scalar ratio are not in complete agreement with the 
Planck 2015 results [27, 28]. 

In this work we showed that using potentials steeper 
than the exponential we can solve the problem of mod¬ 
els of quintessential inflation which give rise to numerical 
values of r larger than the Planck bounds [27, 28], and 
we can obtain a remarkable agreement with the Planck 
2015 results [27, 28]. We have shown that it is possible to 
reproduce the correct post-inflationary evolution during 
radiation and matter eras. The nonminimal coupling be¬ 
tween the scalar field and the neutrinos is shown to drive 
the late-time cosmic acceleration. To summarize, we pre¬ 
sented a successful model of quintessential inflation that 
can describe the entire history of Universe evolution in a 
unified framework. 
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